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A simple and completely general representation of the exact exchange-correlation functional of 
density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any 
Coulomb-interacting system. This representation leads to an alternative point of view on popular 
hybrid functionals, providing a rationale for why they work and how they can be constructed. A sim- 
ilar representation of the exact correlation functional allows to construct fully non-empirical hyper- 
generalized-gradient approximations (HGGAs), radically departing from established paradigms of 
functional construction. Numerical tests of these HGGAs for atomic and molecular correlation 
energies and molecular atomization energies show that even simple HGGAs match or outperform 
state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry. 
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I. INTRODUCTION 

The modern understanding of materials is based on the 
detailed quantitative description of their electronic struc- 
ture afforded by quantum-mechanics. Since electron- 
electron interactions play an essential role in shaping 
material properties, electronic structure must be dealt 
simultaneously with the many-body physics of interact- 
ing electrons. At the heart of this interface between 
many-body theory and electronic-structure theory is the 
exchange-correlation (xc) functional of density-functional 
theory (DFT) [I, 2, s', 4^. 

This functional, Exdn], comprises the exchange en- 
ergy which is known in terms of single-particle 
orbitals, and the correlation energy £'c['t-], which is un- 
known, and must be approximated. An intense and 
multidisciplinary research effort is thus directed at de- 
veloping ever better approximations for i?a;c[n] or -Eci^^]- 
While the local-density approximation (LDA) is still 
widely and successfully used for band-structure calcula- 
tions, total-energy calculations require at least the ac- 
curacy afforded by generalized-gradient approximations 
(GGAs) [6i,iIijiS]- Quantum-chemical applications of DFT 
to atoms and molecules make increasingly use of hybrid 
functionals [l^, [HI mixing a certain (often empirical) 
fraction of exact exchange with LDA or GGA. 

To make progress beyond the GGA level, the con- 
cepts of meta-GGA (MGGA) [H, [H and hyper-GGA 
(HGGA) [a, [11,13 

have been proposed for functionals 
employing the kinetic-energy density or the exchange- 
energy density, respectively, but few explict examples of 
such functionals have been constructed. Established pro- 
cedures of functional construction do not fully determine 
the form of GGAs, MGGAs and HGGAs, enhancing the 
need for introducing empirical parameters. 

Here we propose a hitherto unexplored mode of func- 
tional construction, which provides additional insight 
into the structure of the exact xc functional and the na- 
ture of common approximations to it, and naturally gives 



rise to non-empirical hyper-GGA functionals. 

In Sec.|TT]we use the universal Lieb-Oxford (LO) bound 
[l7| on the exchange-correlation energy of Coulomb- 
interacting systems to derive an exact representation of 
the universal exchange-correlation functional of DFT. As 
an immediate consequence, we also obtain an LO-based 
exact representation of the universal correlation-energy 
functional. 

In Sec. mil the LO-based representation of i?a:cM is 
shown to lead naturally to a reconstruction of the generic 
global hybrid functional, providing additional insight into 
the meaning of the components of the hybrid. Local hy- 
brids are shown to be describable in a similar way. 

Section |TV] is devoted to using the LO-based represen- 
tation of Ec [n] to construct a family of orbital-dependent 
correlation functionals that in the Jacob's ladder clas- 
sification scheme 4] belong on the hyper-generalized- 
gradient approximtion (HGGA) rung. The resulting 
HGGA functional is of different form from other func- 
tionals on the same rung, and satisfies exact constraints 
such as scaling properties and recovery of the gradient 
expansion for weakly varying densities. 



In Sec. IV Al we report numerical tests of the con- 
structed HGGA for atomic and molecular correlation en- 
ergies and molecular atomization energies, and compare 
to common correlation functionals of the LDA, GGA and 
meta-GGA type. We also present, in Sec. IV Bl two fur- 
ther LO-based HGGA correlation functionals, which per- 
form better than that of the preceding section for specific 
properties. One of these employs an unusual multiplica- 
tive self-interaction correction, the other enhances the 
possibilities for error cancellation in combination with 
semilocal approximations for exchange. 

Section IVll contains our conclusions. 
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II. LIEB-OXFORD-BASED 
REPRESENTATIONS OF THE CORRELATION 
AND EXCHANGE-CORRELATION ENERGIES 

The starting point of our analysis of Ec is the 
Lieb-Oxford (LO) bound 1/1], according to which the 
exchange-correlation energy obeys 

E,M>^E^''>], (1) 

where the Lieb-Oxford value for A is Alo = 2.27. Recent 
numerical evidence suggested that the LO bound can be 
further tightened (l8l . [ig| and generalized to one and two- 
dimensional systems [20| . We therefore also report results 
obtained from the value Xel — 1.9555 w 2.0, which is the 
ratio of E^c to E^^^ in the extreme low-density limit of 
the uniform electron liquid (EL), which was conjectured 
to provide the tightest universally applicable bound [19, 

In the form ([T|), the LO bound plays a key role in the 
development of some GGAs [1, [2l|, |23| and meta-GGAs 
[T^ ]. The form of these functionals is dictated by 
other considerations, but typically the value of an oth- 
erwise undetermined parameter in them is chosen such 
that the bound is obeyed for all possible densities. Here 
we make rather different use of the LO bound: instead 
of using it to fix the value of a parameter in a functional 
whose form is obtained in other ways, we use it to deter- 
mine the form of the functional itself. 

By combining the variational principle with the LO 
bound in its form ([1]), we find immediately 

> S,[7i] > E.rAn] > A£;^^^M. (2) 

Exc is thus bounded from above and below, which allows 
us to cast it as 

E^c[n] = (1 - P[n])E,[n] + f3[n]XE^''^[nl (3) 

where P[n] is a density functional taking values in the 
interval [0, 1]. By subtracting Ex[n] we find for the cor- 
relation functional the exact representation 

EM^m{^E!^''>]-EM), (4) 

Since these representations of Exc and Ec are completely 
general, all exact constraints on these functionals become 
constraints on f3[n], with the difference that while ii^ajcM 
and Ec[n\ have values ranging, in principle, from zero to 
—oo, (3[n] ranges from zero to one. 

Although in principle /3[n] is as complex as Exc[n\, we 
note that the maximum absolute error one can make in 
approximating a quantity varying from to — oo is oo, 
while for a quantity varying from to 1 is is 1/2. This 
observation suggests that it may be useful to develop sim- 
ple models for /?[n], designed to recover as many exact 
properties as possible, in order to develop new approxi- 
mations for i^ajcM. 



III. CONNECTION TO GLOBAL AND LOCAL 
HYBRIDS 

Representations ^ and (|3|) have many interesting 
properties, of which we now explore a few. As a first 
application, we compare Eq. ([3]) with hybrid function- 
als, of which a typical single-parameter example can be 
written as 0, [13113 

E'^y" N = (1 - a)Ex [n] + aE^^^ [n] + E^pp^"^ [n] . (5) 

Here a is a (normally empirical) constant determining 
how much LDA exchange is mixed into exact exchange. 

Expression ^ can be considered an approximation to 
the exact representation ([3]), consisting of three steps: 
(i) Replace the functional /3[n] S [0, 1] by a parameter 
a g [0,1]. (ii) Replace the LO parameter A by unity, 
in the second term. According to the general LO bound, 
this means that the correlation energy is underestimated, 
i.e., the resulting energy is not low enough, (iii) This 
underestimate is compensated by adding an explicit cor- 
relation functional £:«pp'-°^ < Q. 

This re-construction of the generic hybrid (O starting 
from the exact representation ([3]) suggests an alternative 
interpretation of the individual contributions to the hy- 
brid functional: E^pp^°^ is not an approximation to the 
full correlation energy, but only to the part missed by 
replacing A = 2.27 by 1. We stress that this is merely 
a change in perspective, as the final form is exactly the 
same. However, such a change may be useful in selecting 
correlation functionals to be used in conjunction with ex- 
change hybrids, and in the construction of novel hybrids. 

A related class of functionals, so-called local hybrids, 
are of the generic form [2^ 

E'J'c^'in] = I dV[(l - a{r))ex[n]{r) + a(r)e^^^[n](r) 

+e^PP™-[n](r)],(6) 

where lower-case letters indicate energy densities. These 
functionals, too, can be interpreted as particular ap- 
proximations to the general LO-based representations, 
by starting from the local LO bound, ea;c[?T.](r) > 
Ae^^"*[7i](r). This local form of the bound (satisfaction 
of which guarantees satisfaction of the global one) is that 
also used in constructing PBE-GGA and TPSS meta- 
GGA. In terms of this bound, the different components 
of the local hybrid can be interpreted in the same way as 
for the global hybrid. However, we note that unlike the 
global LO bound, the local one is not a rigorous prop- 
erty of quantum mechanics, but may be violated. Thus, 
in this sense, local hybrids are less tightly connected to 
universal bounds than global hybrids. 

IV. CONSTRUCTION OF A HYPER GGA 

As a second application, we use representation ^ to 
construct a class of non-empirical hyper-GGAs by en- 
forcing constraints on (3[n\. Occasionally, the expression 
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hyper-GGA is meant to refer to any functional employ- 
ing exact exchange. In this sense, the global and local 
hybrids just discussed are already hyper-GGAs. How- 
ever, we here adopt a more restrictive definition, in which 
hyper-GGA refers specifically to correlation functionals 
that use the exchange-energy density as an ingredient. 

A first approximation to the functional (3[n\ is obtained 
by requiring that the resulting Ec[n\ has the correct uni- 
form density limit. On uniform densities n the exact and 
general representation ^ becomes 



E,[n] = (3[n] {\E^ 



LDA\ 



EM) 



(7) 



and since by definition Ec[n\ — E^ [n] and Ex[n] 



E. 



LDA\ 



we find 



m = 



E 



LDA\ 



{\-l)E^'^^[ny 



(8) 



The use of this /3 in Eq. ^ also for nonuniform densities 
leads to a LO-based functional that correctly recovers the 
uniform gas limit. 

In the same way, we can build in the gradient ex- 
pansion for weakly nonuniform densities n{r), where the 
tilde means that the density is such that the low-order 
gradient-expansion approximation (GEA) is adequate. 
On such densities representation ^ becomes 



E,[fi]^P[fi] (XE^^'^'ifi] - E,[ 



(9) 



and since by definition Ec[n] = E^-^^[h] and Ex[n] = 
E^^^[n], we can construct an approximation to P[n] 
from any functional that on weakly varying densities cor- 
rectly reduces to the gradient expansion. The simplest 
choice would be the gradient expansion itself, leading to 



E^^^[f. 



XE^DA^f,^ - EOEA[f,y 



(10) 



The use of this f3 in Eq. ^ also for strongly nonuniform 
densities leads to an explicit functional correlation func- 
tional recovering the uniform and the weakly nonuniform 
limits, in addition to the LO bound. 

However, the denominator of this prefactor can have 
zeros at some values of the reduced gradient s. The pref- 
actor (3[n\ diverges at such densities, in contradiction to 
the Lieb-Oxford bound. Thus, this bound requires to 
construct f3 from a functional that has the correct gradi- 
ent expansion at small gradients, and does not produce 
divergences at large ones. This problem is solved by using 
GGA instead of GEA, i.e.. 



E^'^^if 



XE^DA^f 



EGGA^f^y 



(11) 



which is also exact for weakly varying densities. The use 
of this P in Eq. ^ for arbitrary densities leads to the 
approximate functional 



e: 



GGAt 



XE^DA[n] ^ EGGA^n^ 



{XEi;^^[n]-EM) 



(12) 



By construction, this functional recovers the uniform 
limit and the gradient expansion to the order built into 
GGA. This is an important caveat, since actually very 
few GGAs recover the GEA both for exchange and cor- 
relation. In fact, the requirement that both the exchange 
and correlation functional used in constructing (3[n\ re- 
produce the low-order gradient expansion excludes pop- 
ular GGAs such as PBE, PBEsol and BLYP, and almost 
uniquely singles out PW91 as the only widely used GGA 
suitable for the construction (2^ . 

Numerically, we have explored many other variations, 
employing different choices of ingredient functionals of 
(3[n\ that do not fully recover the gradient expansion, 
such as PBE. This empirical analysis confirms that the 
choice of PW91 in (3[n] is near-optimal (and certainly 
better than PBE) independently of, but in agreement 
with, the above construction based on recovering exact 
constraints. 

With this choice, Eq. has become an explicit corre- 
lation functional expressed in terms of other known den- 
sity functionals, in particular the exact exchange func- 
tional E^. For this reason, it belongs into the class of 
hyper-GGAs. Interestingly, while representations ^ and 
(|4]) thus rather naturally lead to a connection with hy- 
brids and to hyper-GGA type functionals, they do not 
involve any explicit use of kinetic energy densities, i.e., 
the present functionals belong to the fourth rung of Ja- 
cob's ladder [1, [3, without having passed through 
the third (meta-GGA) rung. 

We note that the LO bound is incorporated in the 
correlation functional through its structure, not by 
choice of a parameter, as in common GGAs and MGGAs. 
It also makes use of the bound for correlation, unlike 
PW91 GGA, PBE GGA and TPSS meta-GGA, which 
use it for exchange. Moreover, Eq. (|12p employs only the 
global LO bound, which has been rigorously proven for 
all possible densities [l3|, and does not make use of the 
local LO bound. 

We also observe that under Levy coordinate scaling 
[13 n(r) ^ n^(r) = -fMl^), both E^^"^ and E^^^^ 
correctly scale as Ej;[nj] — jEx[n]. Thus, on scaling the 
present E^GGAi^^-^ ^j^g exchange scaling factors cancel, 
and E^GGAi j^^j appropriately inherits all scaling proper- 
ties of the correlation functional used in (3. 

Overall, Eq. (fT^ is a nonempirical (in the sense of not 
containing any fitted parameters) correlation functional 
that recovers the uniform and the weakly nonuniform 
limits by construction, correctly scales as a correlation 
functional, and obeys the LO bound for structural rea- 
sons. 

On the down side, in its present form the functional is 
not size consistent. However, we note that size-consistent 
versions of all our functionals can be constructed by re- 
placing the component energy functionals by their re- 
spective energy densities, and integrating over the result, 
very similar to the step leading from a global hybrid to 
a local hybrid. In the specific case of HGGAl, this leads 
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TABLE I: Comparison of our HGGAl functional with standard LDA, GGA and MGGA correlation functionals. First row: 
negative correlation energy of the H atom, in milliHartree. Second row: mean absolute relative error (mare) of the correlation 
energy of 17 atoms, from He {Z = 2) to Ar {Z = 18). Third row: mare of the correlation energy of 35 molecules for which 
highly precise correlation energies are available [S^l- Fourth row: correct for the electron liquid (y/n). 





LDA 


PW91 


PBE 


LYP 


TPSS 


HGGAl(Ai,o) 


HGGAI(Alo) 
MSIC 


HGGAI(Abl) 


HGGAI(Abl) 
MSIC 


H atom (mH) 


21.66 


6.33 


5.71 


[0] 


[0] 


6.24 





6.20 





atoms (mare %) 


119.7 


4.9 


6.8 


3.9 


5.4 


4.75 


4.38 


4.70 


4.37 


molecules (mare %) 


102.5 


7.4 


9.7 


6.7 


9.0 


7.0 


6.6 


6.8 


6.5 


electron liquid 


Y 


Y 


Y 


N 


Y 


Y 


Y 


Y 


Y 



to 

r egg^[n](r) (A4^^[n](r) - e.[n](r)) 
J \eLDA^n]{r) - eGGA^n]ir) ' ^ ' 

which is size consistent (but, as any other approximate 
density functional available, not guaranteed to be size 
extensive) . 

V. NUMERICAL TESTS 

A. HGGAl correlation functional 

In this section we present numerical tests of our 
HGGAl functional. Since this is a correlation functional, 
the most direct and stringent test is the calculation of 
atomic and molecular correlation energies, to which we 
turn first. 

A selfconsistent implementation of orbital function- 
als such as Eq. is possible by means of the opti- 
mized effective potential method and its simplifications 
[2^ . or by the scaled selfconsistency approach 28]. Al- 
ternatively, such functionals can be implemented post- 
selfconsistently on LDA or GGA densities. Below, we 
evaluate all componcnt-functionals of our HGGAs post- 
selfconsistently on PW91 densities and orbitals. 

Table U compares our HGGAl correlation functional to 
a set of widely used correlation functionals: LDA, PW91 
GGA, PBE GGA, LYP GGA, and TPSS MGGA, for 
18 atoms, 35 molecules and the electron liquid. All DFT 
calculations were performed with GAUSSIAN 03 [29,] using 
the aug-cc-pVQZ basis set (except for the molecules from 
Ref. [30|, which uses cc-pVTZ). 

The first row of Table U reports correlation ener- 
gies of the hydrogen atom, which we display sepa- 
rately from those of other atoms because they are ex- 
clusively due to self-interaction, and thus permit to as- 
sess the self-interaction error. We note that already 
HGGAl, without any self-interaction correction, has a 
lower self-correlation error than LDA and PW91 GGA. 
HGGAl-MSIC, which by construction has zero self- 
correlation error, is explained in Sec. IV Bl below. 



The second row reports mean absolute relative errors 
(mare) for atoms He to Ar. As benchmark data we used 
the standard set of CI atomic correlation energies of Ref. 
[13. HGGAl performes better than LDA, PW91, PBE 
and even TPSS MGGA, but loses to LYP. LYP, unhke 
all the other functionals tested here, contains empirical 
parameters fitted to the He atom, which explains its su- 
perior performance when applied to isolated atoms. 

The third row reports molecular correlation energies 
for a set of 35 molecules for which highly precise cor- 
relation energies are available (3^ . Encouragingly, we 
find that on this set HGGAl achieves a lower error 
than all tested nonempirical functionals, including the 
highly sophisticated TPSS MGGA. The fourth row re- 
ports whether the tested functional is correct for the 
uniform electron liquid. All correlation functionals ex- 
cept for LYP (which was not designed to be correct in 
this limit) pass this test. Finally, we note that the per- 
formance of HGGAl is systematically improved by low- 
ering the value of A. This improvement is particularly 
encouraging, as a constraint-based functional should in- 
deed deliver better results when the constraint it is based 
on is sharpened. 

While our functional is, by construction, a correlation 
functional, and thus should be, as a matter of princi- 
ple, compared to other correlation functionals, in prac- 
tice it is clearly important to also test its performance 
for molecular atomization energies. As a test set we em- 
ployed the well established set of 20 molecules that was 
used in the original PBE work Q . Atomization energies 
are calculated from total energies, which in turn require 
chosing an exchange functional in addition to a correla- 
tion functional. To compare like with like, we combine 
all correlation functionals included in Table [J with the 
same exchange functional. Since the spirit of HGGA 261 
is to provide a correlation functional to be combined with 
exact (Pock) exchange, we first calculated all total and at- 
omization energies combining E^pp'^°^ with the exact E^- . 
Results are reported in the first line of Table [Til HGGAl 
performs better than LDA, PW91, PBE and BLYP, while 
HGGAl-MSIC also improves on TPSS. When combined 
with exact exchange, the present HGGA functionals thus 
provide, comparatively, the best atomization energies. 
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TABLE II: Atomization energies of the set of 20 molecules used in Ref. '6]. Values are mean average errors in kcal/mol, and 
were obtained from spin-unrestricted calculations. First row: Using E^pp'^°^ combined with exact (Fock) exchange. Second 
row: Using E^^'^''°^ combining each approximation with the corresponding conventional choice for approximate exchange and 
aU HGGAs with E^^^\ 





LDA 


PW91 


PBE 


BLYP 


TPSS 


HGGAl 


HGGAIMSIC 


HGGA2 


with Ex 


30.8 


29.4 


29.3 


29.1 


26.8 


28.3 


23.7 


33.3 


with SJPP''"" 


31.5 


8.0 


7.8 


4.5 


3.1 


9.1 


13.7 


4.0 



Variants exploiting error cancelation and the 
one-electron limit 



On the other hand, it is well known that error can- 
cellation allows to obtain much better atomization ener- 
gies from local and semi-local correlation functionals by 
combining them with approximate local or semilocal ex- 
change functionals, instead of with exact exchange. Re- 
sults obtained in this way are reported in the second line 
of Table HB where LDA, PW91, PBE and TPSS corre- 
lation were combined with the corresponding exchange 
functionals, LYP correlation was combined with B88 ex- 
change, and all HGGA's were combined with PW91 ex- 
change. 

Evidently, all semilocal functionals benefit hugely from 
this error cancellation. The same applies to HGGAl, but 
the degree of improvement is slightly smaller than for 
GGAs and MGGA. However, we note that our choice 
of P[n] in the construction of HGGAl is not unique. 
The particular combination leading to HGGAl was de- 
termined by the requirement to recover the gradient ex- 
pansion, which turns out to produce rather good stand- 
alone correlation energies. 

In order to provide a fair comparison with function- 
als exploiting error cancellation, we constructed an alter- 
native HGGA that has larger correlation-energy errors 
than HGGAl, and does not recover the gradient expan- 
sion, but combines better with semilocal exchange. This 
functional. 



E. 



HGGA2\ 



E. 



LDAt 



n]iXEi;^^[n]-EM) 



XEP^^^ 



EJn] 



(14) 



is denoted HGGA2 in Table [III and found to give better 
atomization energies than all other tested functionals, ex- 
cept for the very sophisticated TPSS MGGA. It should 
not, however, be used to calculate correlation energies, 
for which it performs much worse than HGGAl. This 
behaviour is similar to that of common hybrids, which 
give excellent atomization energies but much worse cor- 
relation energies. 

While HGGA2 is a variation designed to benefit from 
error cancellation with a semilocal exchange functional, 
we can also build in other desirable features directly on 
top of HGGAl. As an example, we build in a novel type 
of self-interaction correction (SIC). Just as the property 
Ex[n] = E^^^[n] serves as an integrated homogeneity 
indicator, the property i?2;[n*^^)] = — [n^^^], where n'"'^' 



is any one-electron density, is an integrated one-electron 
indicator. We build this indicator into HGGAl by multi- 
plving Eq. (O with the multiphcative SIC (MSIC) factor 



Ex 



MSIC 



EhH 



ELDA]^n]+EH[n] 



(15) 



which was designed to yield zero for one-electron den- 
sities. LYP GGA and TPSS MGGA also achieve this, 
but at a price: LYP erroneously predicts zero correla- 
tion energy for any fully polarized system, while the local 
one-electron indicator used in TPSS MGGA [13] recog- 
nizes only one-electron systems with real orbitals, but 
fails for complex (current-carrying) orbitals. Thus, the 
correct entry for TPSS MGGA, and even more so that 
for LYP GGA, must be interpreted with caution, sig- 
naled in Table H] by [0]. The global one-electron indicator 
Ex[n'''^^] = -EH[n'^'^\ used in our MSIC, does not suffer 
from either problem. 

The resulting HGGAl -MSIC functional is a product of 
three factors, one steming directly from the LO bound, 
one (/3) from the electron liquid, and one {F) from 
the one-electron limit. HGGAl-MSIC spoils the recov- 
ery of the gradient expansion achieved by HGGAl, but 
does correctly recover the one-electron limit. (HGGA 
functionals achieving both properties are currently un- 
der investigation in our group.) As Table [Til shows, 
HGGAl-MSIC combines better with exact exchange than 
any of the other functionals, which is the behaviour ex- 
pected from a HGGA. However, when combined with 
semilocal exchange it is still inferior to functionals ex- 
ploiting error cancellation. 



VI. CONCLUSIONS 

All of the above suggests that the present HGGA func- 
tionals (with or without the MSIC factor) deliver com- 
petitive correlation and atomization energies, matching 
or outperforming those from sophisticated state-of-the- 
art functionals. We stress, however, that in spite of this 
encouraging conclusion we consider the present HGGA 
functionals merely as illustrations of the use of our rep- 
resentations ([3]) and (|31) in the construction and analysis 
of functionals, and not as the final word in this regard. 

Future development of other functionals based on the 
same representations (including novel hybrids) seems 
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promising. The way our HGGAs are constructed from 
the global Lieb-Oxford bound represents a radical depar- 
ture from traditional modes of functional construction in 
quantum chemistry and solid-state physics 0, H, M, B S 



i i, M M M M M E d, mi , which is only beginning 
to be explored. 
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